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Abstract. The continuous geomagnetic field survey holds an important 
potential in future prevention of tsunami damages, and also, it could be used 
in tsunami forecast. In this work, we were able to detected for the first time 
Rayleigh and ionospheric acoustic gravity wave propagation in the Z-component 
of the geomagnetic field due to the Japanese tsunami, 2011 prior to the tsunami 
arrival. The geomagnetic measurements were obtained in the epicentral near 
and far-field. Also, these waves were detected within minutes to few hours 
of the tsunami arrival. For these reasons, these results are very encouraging, 
and confirmed that the geomagnetic field monitoring could play an impor¬ 
tant role in the tsunami warning systems, and also, it could provide addi¬ 
tional information in the induced ionospheric wave propagation models due 
to tsunamis. 
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1. Introduction 

The monitoring of the geomagnetic held has proved to be an auxiliary tool in the study 
of tsunami by several authors [Balasis and Mandea, 2007; Manoj et ai, 2011; Utada et al, 
2011; Klausner et al, 2014a]. The response of the ionosphere to tsunamis and the Rayleigh 
and gravity waves induced by them has been studied broadly over the years [Occhipinti 
et al., 2006, 2008, 2010, 2013; Rolland et al, 2010, 2011; Galvan et al., 2012], and modeled 
[Oeehipinti et al., 2011; Kherani et al., 2012], 

It is well-known that post-seismic acoustic-gravity and Rayleigh waves can be observable 
close to the epicenter (in the near-held, within 500 km with velocity up to = 3km/s) due 
to direct the vertical displacement of the ground induced by the rupture. Indeed, the 
rupture, as a Dirac function, have a broad spectrum of energy including both, acoustic 
and gravity waves as presented on the detailed simulational-observational work by Kherani 
et al. [2012] using TEC and magnetic data. In the far-held, tsunamis induce pure gravity 
waves, and additionally Rayleigh waves induce pure acoustic waves [Occhipinti et al., 
2010]. More details could be found in Occhipinti et al. [2013]. 

The theoretical and observational evidences show that tsunamis can generate 
Rayleigh and acoustic gravity waves (AGWs) in the atmosphere/ionosphere by tsunami- 
atmosphere-ionosphere (TAI) dynamic coupling, and this coupling ehect do not modify 
the main frequencies of these waves. The presence of varieties of wavefronts propagating 
with velocity ranging between Rayleigh to AGW velocity was observed and discussed here 
for the hrst time using ground magnetic data for the Japanese tsunami, 2011 in the near 
and far-held (distances above 500 km). These observations reaffirmed the idea that the 
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magnetogram data can be used as a possible tool for tsunami warnings as it was already 
discussed by Klausner et al. [2014a]. 

2. Dataset 

On the 11th of March at 05 : 46 UT, 2011, a powerful earthquake of magnitude 8.9 
generated a catastrophic tsunami which propagated in the Pacihc Ocean. The epicenter 
was centered at 38.3 °N and Long. 142.4 °E, near to the coast of Japan. To study this 
event, we used the minutely magnetogram data from the Z-component from 9 ground 
magnetic stations which are displayed on Table 1. 

3. Methodology 

In the work of Klausner et al. [2014a], the wavelet analysis was proofed to be an alterna¬ 
tive tool in detection of magnetic helds induced by the tsunami propagation. Nowadays, 
the discrete wavelet transform (DWT) has been used in many different works in geophysics 
[Domingues et al.., 2005; Mendes et al, 2005; Mendes da Costa et al., 2011; Klausner et al., 
2014a, b]. The DWT is based on the multi-scale analysis and local regularities of the sig¬ 
nal. Considering 'ijj as the analyzing wavelet, the wavelet coefficients of these transforms 
are dependent of two parameters: the scale a and the central position of the wavelet 
analyzing function translation b. 

Here, the orthogonal discrete wavelet transform was used [Daubechies, 1992]. The key 
aspect of this transform is that the amplitude of the wavelet coefficients can be associated 
to the local polynomial approximation error which can be dehned by the choose of the 
analyzing wavelet. 
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Using the squared modulus of the wavelet coefficients we can reproduce the discrete 
scalogram related to a discrete scale and a position, where the scale is related to a dyadic 
decomposition in levels j, as a = 2^ and the translation is related to discrete position 
b = 2~% at level j with k G Z. 

In this context, we have that a signal f{t) can be represented by the expansion 

OO OO 

j=—oo k=—oo 

where are the wavelet coefficients computed from the inner-product 

4 = j fit) 'fii‘^^t - k) dt. ( 2 ) 

The wavelet transform in level j -|- 1 is given by 

= 2 ^ ^(m - 2A;) c^, (3) 

m 

where g is a high-pass hlter, is the wavelet coefficient at level j + 1, are the scale 
coefficients at level j, and m E Z [see Domingues et ai, 2005, more details]. 

In this work, we use the Daubechies (db2) wavelet function of order 2, and therefore the 
non-zero low hlter values are h = Also, the sampling rate of 1 

min, as consequence of the time resolution of the magnetic data, gives the pseudo-periods 
of the hrst three levels of 3,6 and 12 minutes, [see Klausner et ai, 2014b, more details], 

4. Results and analysis 

Figure 4 displays, on the right side panel, the magnetic observatories geographic dis¬ 
tribution, and on the left side panel, the Tsunami Travel Times (TTT) map. In both 
graphics, the arrows displays the observatories used in three different wave front of the 
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tsunami propagation direction that we use here to study the tsunamigenic disturbances, 
and these directions are denoted by the numbers 1, 2 and 3. Also for guiding purposes, 
we suggest the use of TTT map to determine the approximated tsunami time arrival for 
each magnetic observatories. 

In the sequence, we will present the analysis of each wave front direction case study. 
The wave front direction number 1 is parallel to the Japanese coast in the west-southwest 
direction, and in this case the magnetic observatories used are HAR, KAK, TTK and 
KNY. Figure 1 shows the travel-time diagram (TTD) of a magnetic induced held trends 
due to the tsunami in the direction 1. As done previous by Kherani et al. [2012], the 
magnetograms of HAR, KAK, TTK and KNY (blue color) have been hltered using a 
Butterworth hlter with 10 to 30 minutes period bandwidth to extract the tsunami-related 
magnetic induced helds from the external ionospheric activity. Also, the tsunami wave 
propagation was simulated using the simulation model developed by Sladen et al. [2007]; 
Sladen and Herbert [2008] (black dot-dashed line). The continuous line represents the 
tsunami wave front propagation with velocity of = 200m/s, the red dashed line is asso¬ 
ciated to Rayleigh waves propagation (= 1.3/cm/s), and the dashed lines are associated 
to acoustic gravity waves propagation = 640m/s. It is well-known that tsunami waves 
can induce long-wavelength acoustic gravity waves (AGWs) with velocity up to = 800m/s 
and short-wavelength gravity waves with velocity of = 200m/s [Kherani et al., 2012], and 
additionally Rayleigh waves with velocity up to = 2.1km/s which can induce pure acous¬ 
tic waves [Occhipinti et al., 2010]. In Figure 1, Rayleigh and acoustic gravity waves fronts 
is observed in the presented TTD. It is possible to notice that depending in the distance 
from the epicenter the acoustic gravity waves can be detected within minutes to hours 
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before the arrival of the tsunami wave. By this reason, we can discriminate the effect (on 
the magnetic data) induced by the ionospheric perturbation, from the effect indnced by 
the secondary magnetic field induced directly by the displacement of the oceanic surface 
(the tsunami). 

In Figure 3, panels from bottom (near-held) to top (far-held) correspond to magnetic 
observatories of HAR, KAK, TTK and KNY, respectively. Each panel, from top to bot¬ 
tom, displays the corresponding magnetogram (Z component), and the square wavelet 
coefficients for the three hrst decomposition levels, (d-^)^ for j = 1,2,3. It is possible to 
notice that Figure 3 shows an increase of wavelet coefficient amplitudes (WCAs) before 
of the tsunami arrival which is owing to the contribution from both geomagnetic distur¬ 
bances and tsunami induced ionospheric disturbances. In particular, following features 
are evident from Figure 3; (a) The wavelet coefficients are amplihed after few minutes 
of the earthquake, and (b) At TTK and KNY, amplihcations are noted few minutes be¬ 
fore the tsunami arrival. These features suggest that the observed disturbances in a few 
minutes after the earthquake may be associated to Rayleigh, and the amplihcations few 
minutes before the tsunami arrival may be associated to the acoustic ionospheric gravity 
propagation, respectively. However, we can only come to this conclusion if these identihed 
tsunamigenic disturbances are observed also in the TTD. This kind of analysis makes pos¬ 
sible to say that these disturbances can be classihed to be from seismic and from tsunami 
origin. 

In a near epicentral held (HAR and KAK), the delay between the tsunami and the 
AGWs is always positive, see Occhipinti et al. [2013] for more details. And in the far-held 
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(TTK and KNY), the AGWs follows the Rayleigh waves arrival, and they can be detected 
even with 2 hours in advance of the tsunami arrival. 

The TTD for the direction 2 is shown in Figure 3 in which the AGWs disturbances are 
identihed following the same strategy as described in the context of Figure 1. The wave 
front number 2 is towards the open sea waters, and in this case the magnetic observatories 
used are OTA, GBI and GUA. It can be again said that such identihcation classihes these 
disturbances to be associated to AGWs wave propagation (= 700m/s) in which appear 
within a few minutes to hours before the tsunami arrival. 

In Figure 5, the vertical component and corresponding three discrete wavelet coefficients 
for the magnetic observatories of OTA, GBI and GUA as indicated in Figure 4 for the 
direction 2. We note from Figure 5 that the magnetic disturbances within a few minutes of 
the earthquake are amplihed at all observatories, these amplihcations are only due to AGW 
propagation. In OTA and GBI, the WGAs around 6 : 00 UT may or not be associated 
to the Rayleigh wave propagation. However, these features related to the Rayleigh wave 
propagation were not detected in the TTD (Figure 3). Also at GBI, the tsunamigenic 
disturbances after the tsunami arrival are associated to AGW propagation. And these 
waves have propagation characteristics similar to the tsunami wavefront propagation. For 
example, the predicted time of the tsunami’s arrival by NOAA at Ghichijima Island was 
at 07 : 14 UT. The tide gauge measurements registered the tsunami maximum height at 
07 : 46 UT with amplitudes up to 1.80m. The GBI magnetic observatory shows remarkable 
wavelet coefficient amplitudes = 1,2,3) due to the tsunami after the tsunami 

arrival time as well. These increase in WGAs correspond to the AGWs propagation. Also, 
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Utada et al. [2011] discussed the similarity between the CBI Z-component fluctuations and 
the IPM fluctuations reported by Manoj et al. [2011] due to the AGW point of view. 

The TTD for the direction 3 is shown in Figure 5 in which the AGWs disturbances 
(= 800m/s) are identified following again the same strategy as mentioned in the context 
of Figure 1. In this case, the wave front number 3 is parallel to the Japanese coast, and 
the magnetic observatories used are ESA and MMB. From Figure 7, it can be again said 
that the magnetic disturbances due to AGWs within a few minutes of the earthquake are 
amplified at all observatories, while the amplification due to AGWs could be see at MMB 
in 2 hours in advance of the tsunami arrival. The reason for that is the ESA observatory 
is near to the epicenter, and an increase of WGA is almost simultaneous to tsunami wave 
arrival. On the other hand, near to the magnetic observatory of MMB, the initial phase 
of the tsunami started at 06 : 38 UT and reached the maximum height (0.74 m) at 06 : 57 
UT. The {d^y{j = 1, 2, 3) showed WGAs around 06 : 00 UT which may be due to Rayleigh 
wave propagation, and again between 07 : 00 and 08 : 00 UT which might be due to AGW 
propagation. At 06 : 57 UT, when tsunami reached the maximum height, the presence 
of an increase in the coefficient amplitudes around this time might be due to tsunami 
maximum height too. 

In all the analysis presented above, post-seismic acoustic-gravity and Rayleigh waves 
were observable close to the epicenter (near-field, within 500 km with velocity up to 
= 2km/s) due to direct the vertical displacement of the ground induced by the rupture. 
Indeed, the rupture, as a Dirac function, have a broad spectrum of energy including both, 
acoustic and gravity waves. Here, the presence of varieties of wavefronts propagating with 
velocity ranging between Rayleigh to gravity wave velocity was observed as presented 


DRAFT 


August 4, 2015, 12:46ani 


DRAFT 



X - 10 


KLAUSNER ET AL.: TSUNAMI EEFECTS ON THE GEOMAGNETIC FIELD 


on the detailed simulational-observational work by Kherani et al. [2012] using TEC and 
magnetic data. However, the Rayleigh wave propagation was only detected here in direc¬ 
tion 1, while AGWs were detected in all directions. This assymetry in the distribution 
of Rayleigh wave propagation was also observed by Galvan et al. [2012]; Kakinami et al. 
[2013] in the TEC data. Their studies also found the Rayleigh waves propagating only in 
the southward direction, although these waves should propagate in all directions. 

5. Conclusions 

Using hltered magnetograms of the geomagnetic Z-component, and also, the spectral¬ 
time diagrams, we presented a detailed TTD with a presence of Rayleigh and AGWs 
which was not done before for the Japanese tsunami using magnetogram data, 2011. 
These perturbations presented velocities of 1.3 km/s and 600 - 800 m/s, respectively, with 
also it is consistent which tsunami-atmosphere-ionosphere coupling theory. These results 
are in conhrmation with the previous studies [Occhipinti et al, 2006, 2008, 2010, 2013; 
Rolland et al., 2010, 2011; Galvan et al, 2012], however they show new insight into the 
tsunamigenic magnetic disturbances due to Rayleigh and AGWs on near and far-held 
distances from the epicenter. 

As explained by Kherani et al. [2012], the presence of periodicity of order of 8 to 30 
minutes in the spectral-time diagrams and, Rayleigh and acoustic gravity waves fronts 
in the TTDs are the aspects of this work which distinguish the magnetic disturbances 
of tsunami in the ionosphere from other external magnetic perturbation sources. Also, 
Kherani et al. [2012] have presented the TTD and reported the presence of acoustic 
disturbances that remained conhned to < 5° which is 500km distance from the epicenter. 
However, they did not see any acoustic disturbance beyond this distance. In this work. 
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we presented these disturbances beyond 500 km from the epicenter distance and up to 3 
hours of the tsunami arrival. 

Also, the fact of Rayleigh waves and AGWs can be detected with a few hours in ad¬ 
vance of the tsunami is very encouraging, and shows that the constant monitoring of the 
geomagnetic field could play an important role to tsunami forecast and tsunami early 
warning systems. 
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Figure 1. TTT map and map of the geographic localization of the magnetic obser¬ 
vatories. TTT map conrtesy of NOAA/NWS/West Coast and Alaska Tsnnami Warning 
Center. 
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Time (hours) 

Figure 2. Z-component dataset filtered with a bandpass with periods between 10 
and 30 minutes. On the vertical axis, it is shown the distance between the earthquake 
epicenter and the magnetic observatory, and on the horizontal axis, the universal time. 
The continuous line represents the tsunami waves propagation at = 200m/s, the red 
dashed line is associated to Rayleigh waves with velocity of = 1.3/cm/s, and the black 
dashed lines are associated to AGWs with velocity of = 640m/s. 
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Figure 3. Magnetograms and the wavelet coefficients for j = 1, 2, 3 for wave front 
case study in the direction 1. 
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Figure 4. Same as Fig. 1 but for the direction 2. 


DRAFT 


August 4, 2015, 12:46am 


DRAFT 




KLAUSNER ET AL.: TSUNAMI EFFECTS ON THE GEOMAGNETIC FIELD 


X - 19 


Figure 5. 
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Same as Fig. 3 but for the direction 2. 
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Figure 6. Same as Fig. 1 but for the direction 3. 
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Figure 7. Same as Fig. 3 but for the direction 3. 


DRAFT 


August 4, 2015, 12:46am 


DRAFT 


























































KLAUSNER ET AL.: TSUNAMI EFFECTS ON THE GEOMAGNETIC FIELD 


X - 21 


Table 1. Network of geomagnetic observatories used in this work. 


lAGA code 

Station 

country 

Geographic coord. 
Lat.(°) Long.(°) 

GET 

Ghichijima 

Japan 

27.10 

142.19 

ESA^ 

Esashi 

Japan 

39.24 

141.36 

GUA 

Guam 

United States of America 

13.59 

144.87 

HARi 

Haramachi 

Japan 

37.62 

140.95 

KAK 

Kakioka 

Japan 

36.23 

140.18 

KNY 

Kanoya 

Japan 

31.42 

130.88 

MMB 

Memambetsu 

Japan 

35.44 

144.19 

OTAi 

Otaki 

Japan 

35.29 

140.23 

TTK^ 

Totsugawa 

Japan 

33.93 

135.80 


^Observatories belonging to Geospatial Information Authority of Japan (GIS). 
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